The non-perturbative QCD potential 
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The potential between heavy quarks can be rigorously defined in a QCD effective fiefd theories framework. I 
discuss the general situation when the coupling constant may not be considered small and provide an explicit 
expression for the potential valid up to order 1/m . Spin dependent and spin independent terms are expressed 
in terms of Wilson loops, completing an ideal journey started twenty years ago with the classical work of Eichten 
and Feinberg. 



1. INTRODUCTION 

The mass of the h and the c quarks is usu- 
ally considered large enough to treat as non- 
relativistic (NR) heavy-quark-antiquark systems 
made up by them (heavy quarkonia: ?/;, T, B^^ 
...). Therefore, these systems are characterized 
by, at least, three separated scales: hard (the 
mass m of the heavy quarks), soft (the momen- 
tum scale IpI ^ mv^ v <^ 1), and ultrasoft (the 
typical kinetic energy E ^^ mv'^ of the heavy 
quark in the bound-state system). This NR pic- 
ture is supported by the success of traditional po- 
tential models. In this case a description of heavy 
quarkonia systems in terms of a NR Schrodinger 
equation is assumed and the interaction potential 
is shaped in order to reproduce as best as possi- 
ble the data (for some reviews see [0~£|). The 
fact that such a description turns out to work 
reasonably well even with simple ansdtze on the 
potential (like, for instance, the Cornell potential) 
suggests that a NR quantum-mechanical descrip- 
tion of heavy quarkonium systems, may, indeed, 
be appropriate and justified from QCD. 

Fig. nl taken from B, shows the static po- 
tential versus the size of different heavy quarko- 
nium systems. It suggests that heavy quarkonium 
ground states may belong to the region of valid- 
ity of perturbative QCD (the region where the 
static potential has a 1/r behaviour). The per- 
turbative QCD potential has been calculated up 
to order a^ (with respect to the energy level) in 
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1^ (the two-loop static potential) and in Ig] (the 
one-loop 1/m potential and the tree-level l/m^ 
potentials). Non-potential effects show up in per- 
turbation theory at the NNNLO and have been 
investigated in 0,§1- Fig- [^ shows, however, that 
most of the quarkonia states lie in a region where 
the inverse of the size of the system is close to 
the scale Aqcd of non-perturbative physics (i.e. 
in the region where the static potential starts to 
rise linearly in r). In this situation the potential 
can no longer be expressed as an expansion in as , 
but, starting from the seminal work of Wilson 
H, it has been related to Wilson loop operators. 
These operators can be eventually calculated on 
the lattice |1Q] or in QCD vacuum models [O. 

Despite the fact that the non-perturbative 
QCD potential has been investigated for more 
than twenty years [[l2|-[l9[ only recently the com- 
plete (in pure gluodynamics) and correct expres- 
sion valid up to order 1/m^ has been calculated 
||2y,|2^. Here, I will report about these last pro- 
gresses. The theoretical framework is NRQCD 
Hi] and pNRQCD [||,0, which are the suitable 
effective field theories for systems made up by two 
heavy quarks. NRQCD is obtained from QCD 
by integrating out the hard scale m. It is char- 
acterized by an ultraviolet cut-off much smaller 
than the mass m and much larger than any other 
scale, in particular much larger than Aqcd- This 
means that the matching from QCD to NRQCD 
can be done perturbatively, as well as within an 
expansion in 1/m |24 2q|. By integrating out 
the momentum scale m,v one is left with the ef- 
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Figure 1. The size of some heavy quarkonia is 
shown with respect to the Cornell potential /O]/. 

fective field theory called pNRQCD, where the 
soft and ultrasoft scales have been disentangled 
and where the connection between NRQCD and 
a NR quantum-mechanical description of the sys- 
tem may be formalized in a systematic way. We 
will assume that the matching between NRQCD 
and pNRQCD can be performed order by order 
in the l/m expansion. 

In |^,|l| the matching of NRQCD to pN- 
RQCD has been performed at order l/ni^ in the 
general situation Aqcd < mv. This has been 
proved to be equivalent to compute the heavy 
quarkonium potential at order 1/m'^. More pre- 
cisely, a pure potential picture emerges at the 
leading order in the ultrasoft expansion under the 
condition that all the gluonic excitations have a 
gap of O(Aqcd)- Higher order effects in the l/m 
expansion as well as extra ultrasoft degrees of 
freedom such as hybrids and pious can be sys- 
tematically included and may eventually affect 
the leading potential picture (like in the pertur- 
bative regime ultrasoft gluons |0]). The obtained 
expression for the potential is also correct at any 
power in ag in the perturbative regime. In the 
following I will introduce the formalism and out- 
line the derivation of the quarkonium potential 
in pNRQCD. I will discuss the results and briefly 
the non-perturbative power counting. 



2. NRQCD 

After integrating out the hard scale m, we ob- 
tain NRQCD |2|]. Neglecting operators involving 
light quark fields |||l, the NRQCD Hamiltonian 
for a quark of mass mi and an antiquark of mass 
7712 up to 0(l/m^) is given by: 
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7J(°) = I d^xi (n'^n" + B''B°) , (2) 

i/^^'°'=-iy'd3xVt(D2+54^V.B)v, (3) 
iy(o,i) =_^(i,o)(^^^.^^2) (4) 
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— d xd^ gfabcGp,^G^^G^^, (5) 

^(0,2)^^(2,0)(^^^.^^2), (6) 

h(1'1) = - /"d'x {dssi'Hxh + dsv^^ (t4,x^ (Tx 

+d,s4^^T'i^X^T'x + d™V'^TVV'X^TVx) , (7) 

where tp is the Pauli spinor field that annihilates 
the fermion and x is the Pauli spinor field that 
creates the antifermion, zD ~ zV + .gA, [D-,E] 
= D-E-E-D and [Dx,E] =DxE-ExD. 
n" = E" + 0(1/771^) is the canonical momentum 
conjugated to A" and the physical states are con- 
strained to satisfy the Gauss law: D • n°|phys) 
= g(7/'^T°V + X^T"x)|pliys}. The coefficients cp, 



cs, d2 and dg can be found in Ref. |24 2 



and dij {i,j = s, v) in p5| for the MS scheme. 

In the static limit the one-quark-one-antiquark 
sector of the Fock space may be spanned by 
|n;xi,X2)(''^ = Vt(xi)xJ(x2)|7i;xi,X2)(°\ where 
|n; xi, X2)'-'^' is a gauge-invariant eigenstate (up 
to a phase) of _ff(°\ as a consequence of the 
Gauss law, with energy En (xi, X2), and Xc(x) = 
ia^X*{^)- |"-;xi,X2)^''^ encodes the gluonic con- 
tent of the state, i.e. it is annihilated by Xc(x) 
and ip{x) (Vx). The normalization condition is 
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(m;xi,X2|n;yi,y2)^''^ = ^nm ]^ (^^^H^i -Yj)- 
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The positions xi and X2 of the quark and an- 
tiquark respectively are good quantum numbers 
for the static solution |n;xi,X2)^°^; n generically 
denotes the remaining quantum numbers, which 
are classified by the irreducible representations of 
the symmetry group Dooh (substituting the par- 
ity generator by CP). We also choose the basis 
such that T|n;xi,X2)(°) = |n;xi,X2)(°\ where T 
is the time-inversion operator. The ground-state 
energy Eq (xi , X2) can be associated to the static 
potential of the heavy quarkonium under some 
circumstances (see Sec. |^). The remaining ener- 
gies Eii (xi,X2), n =/^ 0, arc usually associated to 
the potentials describing heavy hybrids or heavy 
quarkonium (or other heavy hybrids) plus glue- 
balls. They can be computed on the lattice (see, 
for instance, IPTII and [P8|). Translational invari- 
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ance implies that En (xi , X2 ) = En ' (r) , where 
r = xi -X2. 

Beyond the static limit, but still working or- 
der by order in 1/m, the normalized eigenstates, 
|n;xi,X2), and eigenvalues, i?„(xi,X2;pi, P2), of 

the Hamiltonian H satisfy the equations 
r 

iJ|n;xi,X2) = / (i^Xid^a;2|n;x'^,X2} 

2 

xii;„(x;, x'2; p'l, p'2) n '^^'Hx; - X,), (8) 

2 
(m;xi,X2|n;yi,y2) = ^„m J| ^^^^Xj -y^)- (9) 

Note that the positions xi and X2 of the static 
solution still label the states even if the position 
operator does not commute with H beyond the 
static limit. We are interested in the eigenvalues 
En- Eq corresponds to the quantum- mechanical 
Hamiltonian of the heavy quarkonium (in some 
specific situation) . The other energies En for n > 
are related to the quantum-mechanical Hamilto- 
nian of higher gluonic excitations between heavy 
quarks. Expanding Eqs. (||) and (p|) around the 
static solution we get up to 0(1/?ti ) 

2 
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Ef, [yi,y2) - En (Xi,X2J^ 
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Explicit expressions for the energies En , obtained 
from the above formula, can be found in p^,^ . 

3. pNRQCD 

In the static limit, the gap between different 
states at fixed r depends on the dimensionless 
parameter Aqcd^- In general there will be a set 
of states {jT-us} such that EnJ^{r) ^ mv^ for the 
typical r of the actual physical system. We call 
these states ultrasoft. The aim of pNRQCD is 
to describe the behaviour of the ultrasoft states. 
Therefore, in order to obtain pNRQCD all the 
physical degrees of freedom with energies larger 
than rav^ are integrated out from NRQCD. In 
this context one may work order by order in \/m 
(in particular for the kinetic energy), and the 
calculation of the previous section becomes the 
matching calculation between NRQCD and pN- 
RQCD. 

In the perturbative situation, Aqcd'' *C 1, 
{rius} corresponds to a heavy-quark-antiquark 
state, in either a singlet or an octet configuration, 
plus gluons and light fermions, all of them with 
energies of 0{mv^) W. In a non-pcrturbative 
situation, Aqcd?" ^ 1, it is not obvious what 
{tt-us} is. One may think of different possibili- 
ties. In particular, one could consider the sit- 



uation where, because of a mass gap in QCD, 
the energy spUtting between the ground state 
and the first gluonic excitation is larger than 
mv^ , and, because of chiral symmetry breaking 
of QCD, Goldstone bosons (pions/kaons) appear. 
Hence, in this situation, {nus} would be the ul- 
trasoft excitations about the static ground state, 
which we will call the singlet, plus the Goldstone 
bosons. If one switches off the light fermions 
(pure gluodynamics), only the singlet survives 
and pNRQCD reduces to a pure two-particle NR 
quantum-mechanical system. Therefore, the sit- 
uation assumed by all potential models, may be 
now rigorously derived under a specific set of cir- 
cumstances. 

Here, I shall discuss the pure singlet sector with 
no further reference to ultrasoft degrees of free- 
dom. In this situation, pNRQCD only describes 
the ultrasoft excitations about the static ground 
state of NRQCD. In terms of static NRQCD 
eigenstates, this means that only |0;xi,X2)'-°^ is 
kept as an explicit degree of freedom, whereas 
|n; xi, X2)'*'' with n ^ are integrated out. This 
provides the only dynamical degree of freedom 
of the theory. It is described by means of a bi- 
linear colour singlet field, S'(xi,X2,t), which has 
the same quantum numbers and transformation 
properties under symmetries as the static ground 
state of NRQCD in the one- quark-one- ant iquark 
sector. In the above situation, the Lagrangian of 
pNRQCD reads 
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== S'tfiao -/ls(xi,X2,Pl,P2) jS*, (11 



pNRQCD 



where hg is the Hamiltonian of the singlet, pi = 



— iVxi and P2 = 
expansion up to order l/ra^ 

/ls(xi,X2,Pl,P2) 



It has the following 
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The integration of the higher excitations is triv- 
ial using the basis |n;xi,X2) since, in this case, 
they are decoupled from |0;xi,X2). The match- 
ing of NRQCD to pNRQCD simply consists in a 
renaming of things in a way such that pNRQCD 
reproduces the matrix elements of NRQCD for 



the ground state, in particular the energy. This 
fixes the matching condition 



^o(xi,X2,Pl,P2) = /7s(Xi,X2,Pl,P2) 



(13) 



4. HEAVY QUARKONIUM POTENTIAL 

To express the heavy quarkonium potential in 
terms of Wilson loops is quite convenient for lat- 
tice simulations [|0|. We shall use the follow- 
ing notations: (• • .) will stand for the average 
over the Yang-Mills action, Wu for the rect- 
angular static Wilson loop of dimensions r x 
Tw and ((...)) = {...Wa)/{WQ). We define 
{{Oi{ti)02{t2).--On{tn)))c as the connected Wil- 
son loop with Oi(ti), 02(^2), •■• On{tn) opera- 
tors insertions for Tw 1"^ > ti > t2 > ■ ■ ■ > tn > 
—Tw/'^- We also define in a short-hand nota- 
tion lim = lim lim , where T < Tw is the 

T^oo T^ooTw^oo 

time- length appearing in the time integrals. By 
performing first the limit Tw — *■ 00, the averages 
((. . .)) become independent ofTw and thus invari- 
ant under global time translations. 



Using the matching condition (13) and Eq. 
( pl]| ) we get in terms of Wilson loops |2C|] 



V^'Hr)= lim -HWn), 

1 ^00 ± 
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4 lim / dtt((5Ei(i).5Ei(0))), 
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where Ej = E(xj). 

Let us now consider the terms of order l/m^. 
We define 



^(2,0) 



-fe")(r) + T/f/(r)Li.Si 
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where Lj = r x p^. Analogous definitions held 
for l/t"'^), for which, by using invariance under 
charge conjugation plus mi <-^ m2 transforma- 
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tion, we have Vy'{r) = K^,' ^(r), V}r'{r) 
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V^r'ir) ^ndV^X^ir)^VJ^^^^ir;m2 
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mi). 



Using Eqs. ( |l3| ) and (g^ we get, in terms of 



Wilson loops, 



(16) 



X lim / dtt^{gB\{t)gE{m),, (17) 
V^im^r) = -^\im rdt{{[B,,gE^m)c (18) 



lim / dt{{gB,{t)-gBm))c 



„(1)2 



4 T^oojQ 

-^ lim / dh [ dt2 [ dt:i{t2-hf 

X ((.gEi(ti) • 5Ei(i2)5Ei(i3) • gEi(0)))c 
+ i I Vj. lim / dti / dt2 ih - ta)' 



x((.gEl(ti)5Ei(i2)-ffEi(0))), 



x((.gEl(ti)5Ei(i2)-gEi(0))), 



T rt 



--^lim / d<i / dt2{ti~t2) 



Jo 



x(([Di.,5Ei](ii)gEi(i2)-5Ei(0))), 
+ ^ lim / dti2(([Di.,5Ei](t)[Di.,5Ei](0))), 

-^fv;^lim ^^dtt2^(5El(i)[Di.,.gEi](0)))J 
-i lim rdtt3(([Di.,5Ei](i)5Ei(0))),(ViF(°)) 
+ ifv;^lim ^''dtt3^(5El(t)5Ei(0))),(Viy(°))j 
+ i(V2y(2.o)) _ ^ lim Tdtt^ 

x((gEl(i).9Ei(0))),(V;y(°))(V^,y(°)) 
-rffVabc /rf'x lim g((G^,.(x)G^Ja:)G:ja;))), 



--^ir • ^lim^ / dtt{{gBi{t) x ffEi(O))). (19) 

"^ 

For the V^^^^^ potential we define 

^ ^1 , ' (Li • L2 + L2 • Li) + K^i'i)(r) 
2r"^ 

+T/i^5\\r)Li.S2-<^5\'(r)L2.Si 
+Fi^i)(r)Si.S2 + ys(^V^(r)Si2(f), 

where Si2(f) = 3f • cri f • (T2 — cri • cr2. Due 
to invariance under charge conjugation plus 
mi <-> 1712 transformation, we have VJ^ 's i^) — 

V£^^^^(r; rrii ^^ m2). Using Eqs. ( p^ ) and ( p^ ) we 
get, in terms of Wilson loops. 



V''l'^\r) ==if'F lim / dtt^ 
x((gEl(t).9E^2(0)))e, 



T/(M)(^)^,^ £LiL lim / dtf 

^ 2 T^ooJo 

x((gEl(t).9Ei(0)))., 

K^^'^HO = -^(v^v^^^)) 

/"T /"ti pt2 

-i Hm / dii / dt2 / dig (t2 - 13? 

T^^Jo Jo Jo 

x((.gEi(ti) • .gEi(i2).9E2(f3) • .gE2(0))), 
+i( V; lim / dh [ dt2{ti-t2f 

x((5El(ii)5E2(i2)-ffE2(0)))J 



- v; lim / dii / dt2{h-t2Y 

2\ T^OoJg Jg 



x((3E^(ti)gEi(f2)-.gEi(0)))el 
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x((5El(ii).gE2(i2)-5E2(0))>c 

x((5E^(ii)gEi(i2)-.9Ei(0)))c 
4 lim / dh [ dt2{ti-t2f 

x(([Di.,gEi](ti)gE2(t2)-.9E2(0)))e 
4 lim / dti [ dt2{ti-t2f 

x(([D2.,<?E2](ti)gEi(i2)-5Ei(0)))e 
4 lim / dii2^([Di.,gEi](t)[D2.,ffE2](0))), 




; lim / dtt'{{{g-E\{t)[B2.,gE2m)), 

-((.gE^(i)[Di.,5Ei](0))),}j 

i lim / dif^{(([Di.,5Ei](i)ffE^2(0))>c 
x-(([D2.,.gE2](i).gEi(0))),}(ViFW) 

ifv;^lim^^'^dii3|^^^El(t)5E^2(0))>c 
+ ((5E^(i)ffEi(0)))e}(V^,y("))) 



— lim / dtt"^ 

QT^ocJq 

x((5Ei(i)gE^2(0)))c(V;y("))(V^,T/(°)) 
+ idss+d,, lim ((TfT2'')))5(3)(xi-x2), 

V^^'^\r) ^ ~%tr ■ lim / dtt 
x((gBi(i)xgE2(0))), 



(23) 



.1,1)^^^ ^2c^^ 



lim / di((gBi(i) .562(0))) 



-4(4. + d™ lim ((Tfr2°)))'5^'^(xi-X2), (24) 

T\v —^ 00 
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v(i'i)(r) = ^^^^zf'F lim / dt 



(25) 



((gBl(t)5B^2(0))) - ^((5Bi(t) •.gB2(0)))' 



A critical comparison of the above results with 
the previous literature will be done in the conclu- 
sions. Here, I mention that the potentials (111), 
(111) and (p2|) have been calculated for the first 
time in Ref. Ig^l f^^d |£1[ . Since the potential we 
get here is a well defined quantity, derived from 
QCD via a systematic and unambiguous proce- 
dure, and complete up to order 1/rn?, it is not 
affected by the usual ambiguities (ordering, re- 
tardation corrections, etc.), which affect all po- 
tential models and all phenomcnological reduc- 
tions of Bethe-Salpeter kernels Q . For the same 
reason the above result may be relevant for the 
study of the properties of the QCD vacuum in 
the presence of heavy sources. While the lattice 
data for the spin structure of the potential seem 
to suggest the so-called "scalar confinement" , the 
data for the momentum dependent and spin in- 
dependent part of the potential seem to support 
a fiux-tube picture pp9]. It will be interesting 
to see how consistent these pictures are with the 
new momentum and spin independent potentials, 
once lattice data will be available for them. I 
note that some of them are not simply expressed 
by two field insertions on a static Wilson loop. In 
particular an extended object coming from the 
Yang-Mills sector is required (similar extended 
objects would also show up by taking into account 
operators with light quarks). 

5. POWER COUNTING 

The standard power counting of NRQCD, as 
discussed, for instance, in pn], has been tested in 
the perturbative regime. However, even in this 
regime, due to the different dynamical scales still 
involved, the matrix elements of NRQCD do not 
have an unique power counting in v. In the non- 
perturbative regime the problem of the power- 
counting of NRQCD is still open. The above non- 
perturbative formulation of pNRQCD has trans- 
lated this problem to obtaining the power count- 
ing of the different potentials. This is expected 
to be of some advantage. Since the scale mv has 
been integrated out, the power counting of pN- 
RQCD is simpler. Since all the potentials are 



expressed in terms of Wilson loops, there are or 
there will be direct lattice measurements of them. 
Here, let us only say few words about the 
expected behaviour of the potential using argu- 
ments of naturalness on the scale mv, i.e. as- 
suming that the potentials scale with mv. We 
first consider F^"^. In principle, F^"-* counts as 
mv, but, by definition, the kinetic energy counts 
as mu^. Therefore, the virial theorem constrains 



V(°) also to count as 



In the perturbative 



case this extra v suppression comes from the fac- 
tor Qg ~ w in the potential. Using naturalness, 
V^^'^^/m, scales like m,v^. Therefore, it could 
be in principle as large as F^°). This makes a 
lattice calculation of this potential urgent. Per- 
turbatively, due to the factor a^, it is 0{mv'^). 
For what concerns the 1/m^ potentials, the natu- 
ralness argument suggests that they are of order 
m,v'^. However, also here several constraints ap- 
ply. Terms involving VF^*^^ are suppressed by an 
extra factor v, due to the virial theorem. The 
1 dy(o) 



Gromes relations [Ml 



V, 



(2,0) 



V, 



(1,1) 



2r dr ' ^'^ ^_''^^ 

= 0, suppresses by an extra factor v the combina- 
tion y£g ' — Vj^ 'g . Similar constraints also exist 
for the spin- independent potentials ||lj]. Finally, 
it may be important to consider that some of the 
potentials are 0(as) suppressed in the matching 
coefficients inherited from NRQCD. 



6. CONCLUSIONS AND OUTLOOK 

I have reported about a new derivation of the 
QCD potential in an effective field theory frame- 
work. Explicit expressions for the heavy quarko- 
nium potential up to order 1/m?, valid beyond 
perturbation theory, have been written. 

I now compare with previous results in the lit- 
erature. For the spin-dependent potentials we 
find agreement with the Eichten-Feinberg re- 



sults [13 IJ] (once the NRQCD matching coef- 



ficients have been taken into account) except for 
the 1/toiTO2 spin-orbit potential V£ g . Our re- 
sult gives back the well-known tree-level calcu- 
lation, whereas the Eichten-Feinberg expression 
gives half of the expected result. Moreover, our 
perturbative result fulfills the Gromes relation. 
An analysis of this error, present in the orig- 



inal papers and in several others appeared af- 
terwards, is done in [glj. The spin-independent 
potentials have been computed before only by 
Barchielli, Brambilla, Montaldi and Prosperi in 
fl^ ] (the analysis done in [|7|, which appears to 
be inconclusive, has never been published). We 
agree (once the NRQCD matching coefficients 
have been taken into account) with their result for 
the momentum-dependent terms but not for the 
momentum-independent terms, where new con- 
tributions are found. An approach similar to that 
one presented here has been used in pl| in order 
to derive, from the QCD Hamiltonian in Coulomb 
gauge, the spin-dependent part of the potential 
up to 0{\/m?). Our expression (nO) differs from 
that one used in [^, which, in general, would 
give incorrect spin-independent potentials. How- 
ever, if we take our matching coefhcients at tree 
level and neglect the tree-level annihilation con- 
tributions in the equal mass case, we find essential 
agreement for the spin-dependent potentials. 

I conclude, commenting on two possible de- 
velopments of the present work. First, it is 
worthwhile to explore the possibility of express- 
ing the potentials associated with higher gluonic 
excitations in terms of Wilson loop operators as 
done here for the heavy quarkonium ground state. 
The corresponding quantum-mechanical expres- 
sions are known and have been calculated in 
pi] ]. Second, our results are complete at 0(1/to^) 
in the case of pure gluodynamics. If we want 
to incorporate light fermions, the procedure to 
follow is analogous and our results still remain 
valid (considering now matrix elements and Wil- 
son loops with dynamical light fermions incorpo- 
rated) except for new terms appearing in the ener- 
gies at 0{l/m^) due to operators involving light 
fermions. They may be incorporated along the 
same lines as the terms discussed here. 
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1. V. Zakharov (MPI, Munich): At this 
conference a very big perturbative correction was 



reported by Penin at the next order, which has 
not been explored earlier. Does it affect in any 
way your approach or result? 

A. Vairo: The correction mentioned by 
Alexander Penin was first derived in j^. It arises 
from the leading US corrections (i.e. from the 
NLO in the multipole expansion) to the pNRQCD 
Lagrangian in the perturbative regime. 

2. G. Martinelli (Univ. of Rome and Or- 
say): May you comment on the implementation 
of NRQCD in lattice calculations? 

A. Vairo: The implementation of NRQCD on 
the lattice depends crucially on the power count- 
ing in the non-perturbative regime. Once the 
power counting in v is established, the match- 
ing coefficients, which are counted in as, have to 
be included consistently. Due to the breaking of 
Lorentz invariance on the lattice, new terms, with 
respect to those mentioned here, may arise. 
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